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Abstract 

We propose a renormalisable model based on A 4 family symmetry with an SU{5) 
grand unified theory (GUT) which leads to the minimal supersymmetric standard 
model (MSSM) with a two right-handed neutrino seesaw mechanism. Discrete 
Zg X Zg symmetry provides the fermion mass hierarchy in both the quark and 
lepton sectors, while symmetry is broken to identified as usual R-parity. 
Proton decay is highly suppressed by these symmetries. The strong CP problem is 
solved in a similar way to the Nelson-Barr mechanism. We discuss both the A 4 and 
SU{5) symmetry breaking sectors, including doublet-triplet splitting, Higgs mixing 
and the origin of the y term. The model provides an excellent fit (better than 
one sigma) to all quark and lepton (including neutrino) masses and mixing with 
spontaneous CP violation. With the A 4 vacuum alignments, (0,1,1) and (1,3,1), 
the model predicts the entire PMNS mixing matrix with no free parameters, up 
to a relative phase, selected to be 27r/3 from a choice of the nine complex roots of 
unity, which is identified as the leptogenesis phase. The model predicts a normal 
neutrino mass hierarchy with leptonic angles ^43 ~ 8.7°, 042 ~ 34°, ^23 ~ 46° and 
an oscillation phase 6 ^ ~ —87°. 
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1 Introduction 


The Standard Model (SM), although highly successful, leaves many unanswered questions 
in its wake such as: what (if anything) stabilises the recently discovered Higgs boson 
mass? Are the three known gauge forces unified into a simple gauge group which also 
explains charge quantisation? What is the origin of the three families of quarks and 
leptons and their pattern of masses, mixing and CP violation? Why is CP so accurately 
conserved by the strong interactions? The answers to such questions may help resolve 
longstanding cosmological puzzles such as the nature of dark matter and the origin of 
matter-antimatter asymmetry, both of which are unexplained within the SM. 

In this paper we propose a realistic and fairly complete model capable of addressing 
all the above questions unanswered by the SM. The basic ingredients of our model are 
Supersymmetry (SUSY) together with an SU{5) Grand Unihed Theory (GUT), flavoured 
by an ^4 family symmetry (for a review see e.g. [^). The model is minimal in the sense 
that SU{5) is the smallest GUT group and A 4 is the smallest family symmetry group 
that admits triplet representations. Also, below the GUT scale, the model yields the 
minimal supersymmetric standard model (MSSM) supplemented by a minimal two right- 
handed neutrino seesaw mechanism. The model is realistic in the sense that it provides a 
successful (and natural) description of the fermion mass and mixing spectrum, including 
spontaneous CP violation, while resolving the strong CP problem. It is fairly complete 
in the sense that GUT and flavour symmetry breaking are addressed, including doublet- 
triplet splitting, Higgs mixing and the origin of the MSSM fi term. 

The model also allows a WIMP dark matter candidate due to the conserved MSSM 
R-parity, and permits matter-antimatter asymmetry via leptogenesis involving the two 
right-handed neutrinos. We shall show that the leptogenesis phase is equal to the single 
phase appearing in the neutrino mass matrix, providing a direct link between neutrino 
oscillations and matter-antimatter asymmetry, although we shall not discuss cosmological 
aspects any further in this paper. 

We emphasise the predictive nature of the model in the lepton sector, where the entire 
PMNS matrix is predicted without any free parameters, up to a discrete choice of a 
single phase. Large lepton mixing is accounted for by the seesaw mechanism with 
constrained sequential dominance (GSD) [3]-[^). With a diagonal two right-handed neu¬ 
trino mass matrix Mfi, the dominant right-handed neutrino mainly responsible for 
the atmospheric neutrino mass m 3 has couplings to proportional to ( 0 , 1 , 1 ), 

while the subdominant right-handed neutrino giving the solar neutrino mass m 2 has 
couplings to (z/e, Vt) proportional to (1, 3,1). These couplings, corresponding to the so 
called GSD3 scheme [ 6 |[^, originate from A 4 vacuum alignmentj^ The model consequently 
predicts a normal neutrino mass hierarchy, m 3 > m 2 3> mi = 0 . 

As mentioned above, the lepton sector is controlled by a relative phase which is selected 
to be 27r/3, chosen from the nine complex roots of unity arising from spontaneous GP 
violation of a Zg x Zg discrete symmetry, by a mechanism proposed in [^. Such a 


®CSD4 models have been discussed in 



spontaneous CP violating scenario had been proposed previously in order to account for 

We also employ a Zf discrete 
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the smallness of CP violation in the soft SUSY sector 
R-symmetry (as the origin of MSSM R-parity, as in nil) and a missing partner (MP) 
mechanism [T^ for doublet-triplet splitting as recently advocated for flavoured GUTs 


m 


13 . The model predicts very sparse lepton and down-type quark Yukawa matrices, 


with hve texture zeroes, and Yukawa elements involving simple SU{5) Clebsch-Gordan 
(CG) ratios of 4/9 and 9/2 for the hrst and second families, with nir/mb = 1 for the 
third family, all in excellent agreement with their experimental values run up to the GUT 
scale 
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Quark mixing originates predominantly from a non-diagonal and naturally hierarchical 
up-type Yukawa matrix, provided by the broken Zg discrete family symmetry. Quark 
CP violation, however, comes exclusively from a single off-diagonal element in the down 
Yukawa matrix. By contrast, to excellent approximation, all lepton mixing and CP 
violation originates from the neutrino mass matrix, whose structure is also controlled by 
the ^4 family symmetry and the Zg symmetry via the CSD3 type vacuum alignment as 
described above [^. 


Although there have been many attempts in the literature based on A 4 flavoured SU{5) 
SUSY GUTs (for an incomplete list see e.g. [^), we would argue that none are as success¬ 
ful or complete as the present one. For example, many of the previous models predicted 
mixing very close to tri-bimaximal and are by now excluded. Indeed the present model 
is the hrst one based on CSD3 capable of predicting all the lepton mixing parameters 
consistent with current data on lepton mixing (see also HI). The full literature on 
havoured SUSY GUTs, i.e. which involve a (discrete) family symmetry, is quite extensive 
(for an incomplete list see e.g. (^). The goal of all these models is clear: to address the 
questions left unanswered by the SM. It will take some time and (experimental) effort 
to resolve all these models. However the most promising models are those that make 
testable predictions while being theoretically complete and consistent. 


While there are many different chiral superhelds in this model, indeed almost exactly 
a hundred, it is important to note that we are explicitly presenting a renormalisable 
model. Any “non-renormalisable terms” generated below the Planck scale are required 
to have a specihc well dehned realization through multiple renormalisable terms involving 
heavy messenger helds that can be integrated out around the GUT scale. The respective 
effective theory after they are integrated out is actually more predictive than otherwise, 
with a normal neutrino mass hierarchy, a zero lightest neutrino mass, and all lepton 
mixing angles and GP phases predicted. The model presented here is amongst the most 
predictive and complete SUSY GUTs of flavour, consistent with current data. 

The layout of the remainder of the paper is as follows: in Section we describe the 
superhelds directly related to the SM fermions and neutrinos, as well as their Yukawa 
structures as imposed by the GUT and family symmetries when certain A 4 breaking 
vacuum expectation values (VEVs) are applied; we also perform a global ht to the pa¬ 
rameters of the model and present our predictions for the lepton sector. In Section we 
describe the superhelds that are responsible for breaking the family symmetries and how 
the required A 4 breaking VEVs arise. In Section several aspects related to the GUT 
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Table 1: Superfields containing quarks and lep¬ 
tons and symmetry breaking scalars. 


Table 2: Superfield messengers for the quark 
and lepton Yukawa couplings (and other GUT 
breaking couplings discussed in Section 4.2). 


are discussed, particularly how to break SU (5) and Zf down to the MSSM with R-parity 
in a viable way (i.e. addressing doublet-triplet splitting, the origin of the fi term and 
proton decay). We also discuss the resolution to the strong CP problem. In Section]^ we 
discuss the link between leptogenesis and the oscillation phase in this model. Finally in 
Section]^ we summarise our main results and conclude. Appendix [A| summarises the 
conventions used in this paper, in the basis of 
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2 The Yukawa sector of the model 


The model involves an A 4 x SU (5) CP invariant superpotential at the GUT scale, where all 
symmetries, including CP, are spontaneously broken along supersymmetric flat directions. 
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as discussed in Sections and As already noted, it involves a further Zg x Zg discrete 
family symmetry as well as a Zf discrete R-symmetry. The purpose of this section is 
to describe those aspects of the model pertaining to the Yukawa sector, i.e. relevant for 
understanding quark and lepton masses, mixing and CP violation. The flavour sector 
of the model is very important in our approach, since we make a serious attempt to 
understand and, where possible, predict the experimentally observable fermion masses 
and mixing matrices. 

In Table we show the matter superhelds F, Tj that contain the quarks and leptons, 
as well as the right-handed neutrino superhelds Nf and double seesaw superheld T, all 
of which carry unit Zf^ charge. Apart from the A 4 x SU{5) assignments of F (3,5), 
Ti ~ (1,10), Yf ~ (1,1), under Zg they transform as F ~ 0, T* ~ (5, 7, 0), Yf ~ (7, 8). 
Unlike the rest of the quarks and leptons, the right-handed neutrinos are further charged 
under Zg (as are some of the symmetry breaking scalars). 

In Table we also display the six Higgs superhelds, generally denoted H (but also A) 
which serve to break the SU{5) gauge symmetry. The two light MSSM Higgs doublet 
superhelds Hu and Ha will emerge from and a mixture of Yg and by a mechanism 
discussed later. The superheld ^ which breaks Zg is particularly central to this theory, 
as it is responsible for both right-handed neutrino masses and the up-type quark mass 
hierarchy. Finally we have the 6 i superhelds which break Zg and help to control Dirac 
neutrino masses, and nine A 4 breaking triplet havons generally denoted 0 , with various 
vacuum alignments, responsible for large lepton mixing. 


With these assignments, only the top quark gets a mass from a renormalisable Yukawa 
coupling Y 5 T 3 T 3 (which has Zf^ charge 2 as required for an allowed superpotential term). 
All the other quark and lepton Yukawa couplings must arise through higher order terms. 
This provides the basic reason why most of the SM (or strictly MSSM) Yukawa couplings 
appear to be so small. The observed hierarchy of Yukawa couplings between the three 


families will be explained via a discrete Zg version of the Froggatt-Nielsen mechanism 18 


with powers of the low VEV of ^ controlling the hierarchy in the up-type quark sector, 
and also, in part, the smallness of the down quark and electron. 


In order to enhance predictivity we need the messengers listed in Table which is the 
price we pay for having a renormalisable theory at the GUT scale. We denote these 
superhelds either as fermion messengers, Yj, or scalar messengers, Sj, depending on 
whether they carry similar quantum numbers to, respectively, the quarks and leptons 
(with odd Tjf charge) or the symmetry breaking scalars (with even Zf charge). The 
fermion messengers Yj carry similar quantum numbers to down-type quarks and charged 
leptons (and neutrinos). Scalar messengers Sj have quantum numbers similar to Y 5 (the 
superheld that gives the top quark a renormalisable mass term). The S, messengers 
do not get VEVs, which means we need not consider the ehect of diagrams with S, 
superhelds in external legs to the masses of SM fermions. 


The messengers group themselves in pairs of two superhelds with a renormalisable bare 
mass coupling which respects all the symmetries. Their masses are therefore expected 
to be at or around the GUT scale. Although there will be in general distinct masses for 
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different pairs, for simplicity and because they are all expected to be at a similar mass 
scale, we take masses of all such pairs to be M and set it equal to the GUT scale in our 
numerical estimates. We emphasise that the successful predictions of the model in the 
lepton sector (namely predicting the PMNS matrix) is independent of the specihc values 
of these mass parameters. 


2.1 Up quarks 


Apart from the top quark mass, which originates from a lowest order Yukawa coupling, 
the remaining up-type quark Yukawa couplings appear from higher order terms that 
result from combining several renormalisable terms involving Sj messengers and the GUT 
singlet superheld To be precise, the up-type quark Yukawa couplings arise from S, 
messenger tower diagrams shown in Fig. [T] For example, the most suppressed coupling 
arises from the hrst diagram in Fig. [T] Other less suppressed couplings arise from the 
diagrams where at the base one has the respective TiTj, with a shorter tower leading up 
to The least suppressed coupling, the renormalisable operator responsible 

for the top quark mass, is the last diagram in Fig. 


The effective superpotential responsible for the up-type Yukawa couplings is 


lUup — UijH^TiTj 

The resulting symmetric Yukawa matrix for up-type quarks is 


( 2 . 1 ) 


= v - ( — 

I M 



2.18 


where ^ = {^) /M ~ 0.1. The explicit form of Y“ is given in Eq. 
coefficients Uij, which are 0{1) and, by enforcing GP conservation a 


( 2 . 2 ) 

and includes the 
the GUT scale, 

necessarily real. Thus, the hierarchy of the up quark masses as well as the GKM mixing 
angles are given by powers of Due to the structure of this matrix, any phase introduced 
by (^) can be reabsorbed by appropriate redehnition of the three Tj helds, so Y“ does 
not contain a source of GP violation. 


2.2 Down quarks, charged leptons and flavons 

When considering the Yukawa structures of down quarks and charged leptons we must 
inevitably discuss triplet flavonsj^ The assignments of all the flavons under the family 
symmetries appear in Table Indeed, since the three families of F transform as a triplet 
of A 4 (see Table [^, all TiH^F terms require a contraction with at least one A 4 triplet 
flavon to be invariant. 

® As a point of terminology, we refer to as “flavons” any superfields that are GUT singlets transforming 
non-trivially under the family symmetry and that get VEVs. In particular not only A 4 but strictly 
speaking also Zg and Zg are family symmetries, so we also refer to ^ as a “flavon”. 
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Figure 1: Diagrams responsible for the masses and mixings of the up-type quarks. 
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Figure 2: Diagrams responsible for the masses of the down-type quarks and charged leptons. 


From the diagrams shown in Fig. integrating out the fermion messengers X, which 
acquire large masses as a result of either explicit mass terms or GUT scale Higgs VEVs, 
we obtain effective operators of the form 


Hdown — dasTs- 


p , , rji 

-r + 0224 2 " 


M 


M2 


+ di 2 Ti 




{XiY 


(A24) (-H24) 


F, (2.3) 


where dij are 0 {1) couplings. 


doublets inside and ifjg, as discussed in Section |4.2 


The light MSSM doublet Hd is a combination of the 

hence the ^22 term also leads 


to a relevant Yukawa coupling, 
(discussed in Section]^ is 


The alignment of the respective flavon VEVs of 




(4) = !'e 





(2.4) 


such that, apart from di 2 , the contraction appearing with Ti^ 2,3 isolates the respective 
Fi^ 2,3 family. This would lead to diagonal Yukawa structures if not for the additional 
term connecting Ti(0^F) (see Fig. 2d). 


6 































The resulting effective Yukawa matrices are, schematically: 


/ (0 Ve {0 Vf. 


yiR- 
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V 
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A 24 

0 

0 


^A24^i?24 

M 2 

0 


(2.5) 


M/ 


where UA 24 and VH 24 are the respective VEVs of A 24 and M 24 (given in Eq. 2.6), and we 
include the subscripts LR to emphasise the role of the off-diagonal term to left-handed 
mixing from Y'^. The off-diagonal term in Y® also provides a tiny contribution to left- 
handed charged lepton mixing 0^2 ~ which may safely be neglected. It also 

introduces CP violation to the CKM matrix via the phase of 


Furthermore, because the underlying renormalisable theory is known, the diagrams in 
Fig. [^are the only contributions for each family. The SU{5) contractions and associated 
CG coefficients appearing for each family are unique [I^[l4] . With the GUT scale sym¬ 
metry breaking as discussed in Section each of the scalars here get a VEV with the 
group structure: 


(H-.r = 5 ? v,/V2 

{«j5)f = (5? - 4 “'a" - 44“'4)4' 

{H 2 A)t = cliag(2,2, 2, -3, -3) 
{A24)b = diag(2,2, 2, -3, -3) va^^, 


where the indices run a,b,c = 1,..., 5. This leads to the GUT scale prediction: 
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(2.7) 


The explicit forms of Y'^ and Y®, including GG 
Eq. |2.19| and Eq. |2.20[ respectively. 


and dij coefficients, are given later in 


2.3 Neutrinos and CSD3 

In order to obtain the GSD3 vacuum alignment in this model we couple the neutrinos to 
a set of flavons distinguished by the Zg symmetry. Of the superhelds in Table [T| only 
the right-handed neutrinos and some of the flavons are charged under this symmetry. 
For clarity, we relabel two of the flavon helds as 0atm = 03 and 0soi = 04 , to highlight 
their role in producing neutrino mixing. We also write = Yf to denote the right- 
handed neutrino that dominantly leads to the atmospheric neutrino mass, and = Y| 
as that which contributes mainly to the solar neutrino mass. The relevant terms in the 
superpotential giving neutrino masses are thus 

W'. = SiTfsf (2.8) 
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The flavons 0atm and 0atm gain VEVs according to the CSD3 alignment: 


(0atm) 





(0sol) — "i^sol 



(2.9) 


where Watm and Vgoi are generally complex. Denoting the phases of VEVs as pi = arg(nj), 
only the relative phase Patm —Psoi between the VEVs is physically relevant, and is con¬ 
strained to a discrete set of valnes, as discnssed in Section 3.3 The flavon ^ (already 
responsible for the np qnark masses) is also acting as a Majoron by generating hierar¬ 
chical right-handed nentrino masses. At the effective level, the Dirac terms resnlt from 
conpling the nentrinos (and H^) to 0 atm and 0 soi via the flavon 6*2 (an ^4 singlet carrying 
Zq charge). The corresponding diagrams with associated messengers appear in Fig[^ 
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Figure 3: Renormalisable diagrams leading to neutrino effective terms. Diagrams (a) and (b) 
are responsible for neutrino Yukawa terms (leading to Dirac masses) while (c) and (d) give 
right-handed neutrino Majorana mass terms. 


In tnrn, the Majorana mass term for is also non-renormalisable and we refer to the 
snperheld T as the respective messenger. It conples only to and simply provides the 
non-renormalisable mass term for snppressed relative to the mass of As T has 

the qnantnm nnmbers of a third right-handed nentrino, one can also consider this held as 
mediating a donble seesaw mechanism, responsible for the mass. The mixing term 
]^"^atm-^soi! thongh allowed by the symmetries, is absent as there is no combination of 
messengers able to prodnce it. 


We write {^) = where is chosen from a discrete set of available phases, discnssed 

in Section 4.1 (see Eq. |4.2 ). This phase originates from the spontaneons CP violation of 
a discrete Abelian symmetry l^f^, in onr case the Zg. We will now show that p^ and 


Patm“psoi fix the relative phases within the effective nentrino mass matrix and conseqnently 
the leptonic mixing angles. 


In a Snpersymmetric (SUSY) model, the relevant terms in the snperpotential giving 
nentrino masses, in the diagonal charged lepton basis, are 

w. = + A4..„K, „/V'„ + (2.10) 
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where Lj are three families of lepton doublets and the (CP conjugated) right-handed 
neutrinos A^ftm with real positive masses Matm and Msoi do not mix. 


The structure of X'^ is determined by the vacuum alignments of 0atm and (psoi- The Dirac 
and Majorana matrices, derived by comparing the superpotential terms in Eqs. 2^ and 
I2.ini are 


= 


where a = {O 2 ) and b = SeKsoi/ (flr)- 


/O b\ 


( 2/3 (0^ 

\ 

0 

(a 36 1 

M^ = 

M 

\a b ] 


\ 0 

2/4 (0 / 


( 2 . 11 ) 


For the see-saw mechanism we shall introduce a different convention for Yukawa and 
Majorana masses. The Yukawa matrices Y®, are dehned in a LR convention by 

( 2 . 12 ) 

where i,j = 1,2,3 label the three families of lepton doublets Lj, right-handed charged 
leptons e):j and right-handed neutrinos below the GUT scale; H‘^ are the Higgs 
doublets which develop VEVs Vu, Vd- The physical effective neutrino Majorana mass 
matrix mX is determined by the seesaw mechanism, 

nf = (2.13) 

where the light Majorana neutrino mass matrix rrY is de£nec0by = —^mXV luI+Ii.c., 
while the heavy right-handed Majorana neutrino mass matrix Mr is dehned by = 
-IMrIYrPr + h.c.. 


There is a simple dictionary between the seesaw basis and the SUSY basis, as follows: 
compared to the SUSY basis in Eq. |2.8| used in leptogenesis calculations we see that 
while Mr = {M^y. Hence the neutrino matrices become, in the seesaw 


= (A'' 
basis. 


= (A^)* = 



Mr = {My = 



\ 


ivAorJ 


(2.14) 


Seesaw produces the effective neutrino mass matrix 

/O 0 0\ /I 3 1\ 

m'^ = ma\0 1 1 + 3 9 3 , (2.15) 

Vo 1 i; vi 3 1/ 

where ma = v‘^\a\‘^/{y^lvy/M) and nib = /{y^vy and we have multiplied through¬ 

out by an overall phase which we subsequently drop, keeping only the (physical) relative 
phase 

h — T 2(pa,tm Psol)) (2T6) 


^The conventions for and m'^ differ, respectively, by overall Hermitian and overall complex 
conjugation compared to those used in the Mixing Parameter Tools package [19| , which was used when 
performing global fits. 
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where we recall the above dehnitions of phases, 


Pi = arg((0), Patm-Psoi = arg(watm<oi)> (2-17) 


and that CP conservation at high energies ensnres that Ui and M are real. By argnments 
given in Section |3^ and Section |4T, we can restrict the physical phase rj to a discrete 
choice, namely one of the nine complex roots of nnity. The valnes r] = ±27r/3 are preferred 
by CSD3 [^[^. Note that the model predicts a normal nentrino mass hierarchy, namely 
m 3 > m 2 3> mi = 0 , which will be tested in the near fntnre. 


The sign of r] has phenomenological significance, as it fixes the leptonic Dirac phase 6K 
Specihcally, a positive rj nniqnely leads to negative and vice versa [^. As experimental 
data hints at 5^ ~ l‘^i fhe a posteriori preferred solntion has positive rj = +27r/3. The 
sign of r] also has cosmological signihcance, as discnssed in Section For example a 
positive r] = +27r/3, together with the reqnirement that baryon asymmetry is positive, 
implies that the lightest right-handed nentrino shonld be A^i = A^atnn while 
shonld be somewhat heavier, which is the natnral ordering in onr model. 


2.4 Full parameter fit 


The strnctnre of the Yukawa matrices and neutrino mass matrix is set by the theory, 
up to 0{1) coefficients. The VEVs of the fields A 24 and 7^24 are at or near the GUT 
scale, but otherwise undetermined. This freedom coincides with the choice of coefficients 
in the Yukawa matrices, providing no extra degrees of freedom in the determination of 
the Yukawas other than to provide the appropriate scale. The same is true for the flavon 
fields (pe, (pfj. and cpr, which provide the necessary hierarchy in the down-quark and charged 
lepton Yukawa sector. 


The neutrino matrix m^ is given in Eq. 2.15 Letting Vf represent the VEV of a field /, 
the Yukawa matrices are as follows: 


Y'^ = 






«12 11=^1 

M 22 |f^| 

M23|f| 

(2.18) 


M23|f| 

■*^33 / 



y'd _ _ 


1 


/I, \ViVe\ 

4 


0 

0 


il2 


2 d 


^^24^/^ I 




22 ' 


M 2 




0 

0 


. Yt\ 


(2.19) 
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0 


0 


F" = — 


1 


/ 1 \v^Ve\ 

9 ^Va24P 


\ 


V 


\VA2ivH2. 

0 


i22- 


M2 


0 


0 


c?33^r^ 


w, 

W'J 


( 2 . 20 ) 


As already remarked, the phases in Y'^ from powers of (^) = can be removed 

by held redehnition. Without loss of generality we have rephased helds such that the 
only phase appearing in and is the phase ( as shown in Eqs. 2.19, 2.20, so all 


quark CP violation originates from the single phase ( appearing in Yf 2 - In turn, ( is 
determined by a combination of phases coming from various held VEVs; more precisely, 
C = p^— 2 pH 24 —PA 24 - As long as it is reasonably far from zero, it can produce the necessary 
CP violation. Diherent choices of ( do not ahect the goodness-of-ht, corresponding simply 
to diherent but equally valid choices of 0{1) coefficients. For our ht we choose C, = n/3. 
Note that the corresponding phase in Yg'i does not contribute to leptonic CP violation, 
since this term does not ahect left-handed mixing, to an accuracy of 0 {me/m^). 


To ht the real coefficients dij, m, 


and mb, we propose a function that relates the 
N physical predictions Pi{{x}) for a given set of input parameters {x} to their current 
best-ht values pi and their associated la errors, denoted ai, by 


^ 2 ^^ r mx})-p4 ^ ^ 
i=i h 


( 2 . 21 ) 


The errors Uj are equivalent to the standard deviation of the experimental hts to a Gaus¬ 
sian distribution. For most parameters, this is essentially the case, with the exception of 
the (lepton) atmospheric angle 6 ^ 23 . For a normal hierarchy, the distribution is roughly 
centered on 6*23 = 45°, while the best ht value is given by 6*23 = 42.3°. So as to not 


overstate the error for 623 , we approximate its distribution by a Gaussian about 42.3°, 
setting agi^^ = 1 . 6 °. 

We now wish to minimise In this analysis, = 18, corresponding to six mixing angles 
6 [j (neutrinos) and (quarks), the GKM phase 5'^, nine Yukawa eigenvalues for the 
quarks and charged leptons, and two neutrino mass-squared diherences Am^i and Am^^. 
In the lepton sector, we use the PDG parametrisation of the PMNS matrix 20 Gpmns = 
-^ 23 ^i 3 -^i 2 ApDG in terms of Sij = sin 6 *)-, = cos 6 \-, the Dirac GP violating phase 5^ and 

further Majorana phases contained in Ppoc = diag(l, e*^). Experimentally, the 
leptonic phase 6 ’' is poorly constrained at Icr (and completely unconstrained at 3 (t), so is 
not ht, and left as a pure prediction of the model, as are the (completely unconstrained) 
Majorana phases 021 and 031 . As this model predicts only two massive left-handed 
neutrinos, i.e. mi = 0, one Majorana phase is zero, which we take to be 031 = 0. 

The running of best-ht and error values to the GUT scale are generally dependent on 
SUSY parameters, notably tan/3, as well as contributions from SUSY threshold correc¬ 
tions. We extract the GUT scale GKM parameters and all Yukawa couplings (with associ¬ 
ated errors) from 21 for judicious choices of tan /3. In further reference to 1^, we choose 
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Parameter 

(from [^) 

tan/3 

= 5 

tan/3 

= 10 

IJ-i 

CTi 

Idi 

(Ji 

to 

13.027° 

0.0814° 

13.027° 

0.0814° 


0.1802° 

0.0281° 

0.1802° 

0.0281° 

Oh 

2.054° 

0.384° 

2.054° 

0.384° 

51 

69.21° 

6.19° 

69.21° 

6.19° 

Vu 

2.92 xlO -6 

1.81 xio-6 

2.88 xio-® 

1.79 xio-6 

Vc 

1.43 xlO-3 

1.00 xio-"‘ 

1.41 xio-3 

9.87 xio-® 

Vt 

5.34 xio-i 

3.41 xio-2 

5.20 xio-i 

3.15 xio-2 

Vd 

4.81 xio-6 

1.06 xio-® 

4.84 xlO-® 

1.07 xio-6 

Vs 

9.52 xio-5 

1.03 xio-5 

9.59 xio-^ 

1.04 xiO”® 

Vb 

6.95 xio-3 

1.75 xio-4 

7.01 xio-3 

1.78 xio-4 

Ve 

1.97 xio-6 

2.36 xio-« 

1.98 xio-® 

2.38 xio-8 


4.16 xiO”"^ 

4.97 xio-6 

4.19 xiO”'* 

5.02 xio-6 

Vr 

7.07 xio-3 

7.27 xio-5 

7.15 xio-3 

7.42 xlO”® 


Table 3: Best fit values for quark and charged lepton parameters when run to the GUT 
scale as calculated in 21 , with the SUSY breaking scale set at 1 TeV. We have included an 


overall contribution from threshold corrections corresponding to fji, = —0.24375 which affects 
primarily the b quark Yukawa coupling yi,. y,i represents the best-fit value and Ui the error, as 
defined in Eq. 


2.21 


for the parameter rji, parametrising the threshold corrections a valne f]b = —0.24375; a 
non-zero valne is reqnired primarily to prodnce a necessary (small) difference in b and r 
Ynkawa conplings. Experimental nentrino parameters are extracted from 
is reprodnced in Tables and 


22 . All data 


Parameter 

Mi ± lo- 

3(7 range 

(from [22|) 

oh n 

33.48 

31.29 ^ 35.91 

oh n 

42.3 1?;^ 

38.2 ^ 53.3 

oh n 

Q r +0.20 
0.0 _0.21 

7.85 ^ 9.10 

5^ (°) 

306 

0 ^ 360 

Amii 

10-5 

7.50 

7.02 ^ 8.09 

Amii 

10-3 

+2.457 

+2.317 ^ +2.607 


Table 4: Table of current best fits to experimental data for lepton mixing angles and neutrino 
masses case of normal mass squared ordering taken from 
ranges. 


22 , with 1(7 and 3a uncertainty 


Minimisation by differential evolntion was performed in Mathematica, prodncing the set 
of 0(1) inpnt coefficients and the corresponding physical parameters seen in Table 
with an associated = 7.98 (for tan/? = 5) and = 7.84 (for tan/? = 10). 

In this £t, the VEVs of A 24 , 7^24 and the three 4>e,ti,T are hxed by hand in terms of 
the scale M, which is taken to be the GUT scale, i.e. M ^ 3 x 10^® GeV. Similarly, the 
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Higgs doublet VEV enters only implicitly through nia and nif,, but is understood to take 
(at the GUT scale) the value vh = 174 GeV. We set 


ug = 6 X 10 
= M 

vh24 = 3 X 10“^M 


Ve = 10“^M 

= lO-^M 

= 5 X lO'^M 


( 2 . 22 ) 


The value of is chosen to accommodate not only the £t to parameters but also to 
control the /r-term, as discussed in Section 4.2 Meanwhile the factor ~3 split between 
nA 24 and VH 24 , assists in establishing a hierarchy between the e and /i families. With the 
above numerical values for the VEVs, we get the following Yukawa matrices in terms only 
of 0 { 1 ) coefhcients and the complex phase C,\ 


Y'^ 




ye 


1 

71 

1 

7! 


1.296x10“® ■ uii 

2.16x10“'^ ■ Mi 2 

3.6x10“® ■ Mi 3 

2.16xio-'‘ ■ U21 

3.6x10“® ■ U22 

6x10“® ■ U23 

3.6x10“® • U31 

6x10“® • Us2 

U33 

1.5x10“® ■ dll 

2x10“^^ ■ di2e*^ 

0 \ 

0 

6x10“'* ■ d22 

0 

0 

0 

5xl0“® • d33/ 

6.67x10“® ■ dll 

0 

0 \ 

2xi0“4 ■ di2e*‘’ 

2.7x10“® ■ d 22 

0 

0 

0 

5xl0“® • d 33 / 


(2.23) 


(2.24) 


(2.25) 


It is worth reiterating that the neutrino mass matrix phase rj can be forced to admit only 
phases coming from the nine complex roots of unity, essentially due to spontaneous GP 
violation with the Zg symmetry, where we select r] = ±27r/3, both of which yield equally 
good fils, with only minor adjustments to 0{1) coefficients. The primary effect is in 
the prediction of 5*; as previously stated, positive r] corresponds to negative 6 K As this 
is preferred by experiment, the results presented in Table are for r] = +27r/3. The £t 
also predicts the Majorana phases 021 = 72° and 031 = 0. 


In order to understand the signihcance of the fit; and assess the strength of the model 
overall, it is prudent to enumerate the parameters and predictions of the model. The 
nominal parameter count at the GUT scale is very large, owing to the diverse held content. 
However, at the scale where we are able to make predictions, many of these parameters 
combine to give a constrained set of free parameters that need to be determined. Notably, 
the VEVs of Higgs and havon helds such as those given in Eq. |2.22 do not constitute true 
degrees of freedom, as they can be absorbed by redefining other parameters. 


Relevant parameters that require consideration include: six four dij, masses rua 
and mb, phases rj and (, the threshold factor f]b, and tan/3, for a total of Nj = 16 
input parameters. However three of these parameters, namely tan /3, rj and (, are hxed 
prior to the £t, with the latter two phases restricted to discrete choices, as discussed 
previously. Finally, the factor f]b affects only the coupling yb and is htted by hand. As 
mentioned earlier, the model hts N = 18 observables, including nine Yukawa eigenvalues. 
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Table 5: Fitted input quark Yukawa coefRcients Uij and dij (arranged by their position in the 
Y“ and Y'^ matrices, respectively) and neutrino mass parameters nia and mi), and associated 
physical parameters produced by minimising the function We choose = tt/S. With 
tan/3 = 5, the fit gives = 7.98, while tan/3 = 10 gives y^ = 7.84, both very good fits. The 
largest single contribution to y^ is from the fit to the atmospheric angle 6*23 ■ These results are 
given with rj = +27r/3. The non-zero Majorana phase is also predicted to be 021 = 72°, and 
is insensitive to tan /3, as indeed are all the mixing angles and phases. 


two neutrino mass squared differences, six mixing angles and the quark CP phase. In 
addition the model predicts the leptonic CP phase 6 \ two Majorana phases (one of which 
is zero) and a massless physical neutrino. 


3 A 4 symmetry breaking and the flavon vacuum 


In order to address ^4 symmetry breaking we need to address three aspects of the flavon 
vacuum: what drives some flavons to have VEVs a few orders of magnitude below the 
GUT scale, what determines their vacuum alignment, and what Exes the relative vacuum 
phase Patm —Psoi = ai’g('yatm'ysoi) (a^d consequently the physical phase p). In this section 
we consider each of these issues in turn. 
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3.1 Driving the flavon vacuum expectation values 


The flavon cf) VEVs are driven by radiative breaking (see e.g. 24 for a recent review). 
The soft squared mass terms appear as /Aj) and become negative at the 

required scales Aj < M, lifting the respective flat direction to a few orders of magnitude 
below the GUT scale for each flavon: ~ A^. Hierarchies of VEVs are thus natu¬ 

rally expected due to the logarithmic nature of this mechanism. The precise symmetry 
breaking scale Aj for each field 0 , depends on otherwise undetermined parameters in the 
model which are different for each flavon, such as superpotential terms involving messen¬ 
gers]^ Therefore a hierarchy for such flavon VEVs is generated and remains stable due 
to radiative breaking 


25 . On the other hand those for n, 


g ^ V ^ Vq 


in Eq. 


2.22 


arise from 


F-terms as discussed later. 


3.2 Flavou vacuum aligumeut 


Thus far we have assumed that the A 4 triplet VEVs are aligned in special directions. In 
this section we describe how these directions are obtained by the superpotential terms 
allowed by the symmetries. In doing this, the role of Zg becomes clearer. The driving 
sector, a set of superflelds A* and Oij with Zf charge 2, is listed fully in Table [7|p| The 
inclusion of the Zg symmetry is necessary because the driving superpotential, responsible 
for aligning the flavons, needs to have each driving field isolated, as shown below (see 
also 1 ^): 


hhalign -|- ^2(0202 + 026 ^ 1 ) 

T Geji0e0/j T GeT0e0T T G^t-0^0t- 

+ Ge30e03 + ^230203 + Gi20102 + Oi30i03 (3.1) 

+ + ^25 0 2 05 + 0^60/^06 + O5g050g 

-|- Og 40 g 04 -|- ( 9 i 40 i 04 . 


The additional Zg charges ensure each term is separated from all others, leading to an 
array of vanishing F-term conditions that force mutual orthogonality conditions between 
many of the vacuum alignments. Since this was fully discussed in [^ , we need only state 

^Examples of such a Yukawa coupling for the neutrino flavons are 0atmt^2fi4 and 0soF2fi2. 

®Note that the O (and P) fields that carry no Zg charge couple to (with some power of 

e.g. P 22 H 5 H^. We do not discuss these further as the respective F-terms do not affect the alignment 
nor the origin of the /i-term. 
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Field 

Representation 

Ai 

SU(5) 

Zg 

Zg 



1" 

24 

0 

0 

2 

^2 

1" 

24 

3 

0 

2 

^3 

T 

24 

3 

0 

2 

Ti 

T 

24 

7 

0 

0 

T2 

1" 

24 

2 

0 

2 

T3 

T 

24 

5 

0 

0 

T4 

1" 

24 

4 

0 

2 

Ts 

T 

24 

4 

0 

0 

Te 

1" 

24 

5 

0 

2 

Tj 

T 

24 

2 

0 

0 

Ts 

1" 

24 

7 

0 

2 

Tg 

1 

75 

0 

0 

0 

Tio 

1 

75 

0 

0 

2 

Til 

1 

75 

6 

0 

0 

Ti2 

1 

75 

3 

0 

2 


1 

50 

4 

0 

2 

D 2 

1 

50 

3 

0 

0 


1 

50 

1 

0 

2 


1 

50 

8 

0 

0 

Hi 

1 

75 

6 

0 

2 

TI 2 

1 

75 

3 

0 

0 


Table 6: Superfields that break R symmetry 
and cause doublet-triplet splitting through the 
missing partner mechanism. 


Field 

Representation 

d4 

SU(5) 

Zg 

Zg 


A, 

3 

1 

3 

0 

2 

Ar 

3 

1 

4 

0 

2 

A2 

3 

1 

7 

0 

2 

Oefi 

1 

1 

6 

0 

2 

Oer 

1 

1 

2 

0 

2 

O^r 

1 

1 

8 

0 

2 

Oe3 

1 

1 

6 

5 

2 

O23 

1 

1 

5 

2 

2 

O12 

1 

1 

5 

1 

2 

Ol3 

1 

1 

3 

3 

2 

Ofj.5 

1 

1 

0 

4 

2 

O25 

1 

1 

2 

1 

2 

O^G 

1 

1 

1 

4 

2 

O5G 

1 

1 

7 

2 

2 

Og 4 

1 

1 

2 

3 

2 

Ol4 

1 

1 

4 

3 

2 

Pee 

1 

1 

0 

0 

2 

PtJ.il 

1 

1 

3 

0 

2 

P22 

1 

1 

7 

0 

2 

Pe4 

1 

1 

7 

5 

2 

Pie 

1 

1 

6 

4 

2 

P44 

1 

1 

5 

4 

2 

T34 

1 

1 

4 

4 

2 

pl,2 

-^33 

1 

1 

3 

4 

2 

P2r 

1 

1 

1 

3 

2 


Table 7: Driving superfields for the flavon 
alignment potential. 


the results here, namely that Eq. 3.1 leads to 


the following vacuum alignment patterns: 















(3.2) 


The role of the VEVs (containing two zero entries) of the superhelds (pe,fi,T appearing in 
Eq.|2 ^ was already discussed in Section |2.2[ Meanwhile, the role of the VEVs of the 
flavons 03,4 (redubbed 0atm,soi) was described in Section 2.3, It is the special structure 


of these vacuum alignments, combined with the phase of rj in the neutrino mass matrix. 
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that leads to the very successful prediction of the leptonic mixing angles (as described in 
Section 2.4). The remaining VEVs are not directly relevant to the masses and mixings 
of SM fermions, but help shape the VEVs of 0atm and 0soi- 


3.3 Flavor! vacuum phases 


With the direction of the A 4 triplet flavons (j) hxed, we turn now to a discussion of 
how to £x the relative phase Patm —Psoi = arg(uatmnsoi) to a discrete choice. We present 
a mechanism which does this by adding a number of fields P that are A 4 and SU{5) 
singlets, also given in Table and which resemble the O helds except they do not force 
orthogonality between the flavons 0. 


These fields and their respective charge assignments result in the following invariant 
superpotential terms: 

Wphase ~ -Pee(0e0e + + P^e) + PlifJ.i.^’tJ.4’ + ^2^3 + P^^) 

+ -Pe4(0e04 + 6^102) + -^22(0202 + 

+ -Ple(010e + 060 t) + ^44(0404 + 050r) + ^34(0304 + 060e) 

+ -^33 (0303 + 010/i + 050e) + -P2r(020r + 0306 + 0405); 

where each term technically has an associated real coupling A which is 0{1) and may be 
made positive by held redehnitions. We omit these for simplicity as they have no effect 
on the general argument presented here, with one caveat: the two superhelds have 
exactly the same quantum numbers but different A couplings to havons. Due to this du¬ 
plication there are two independent relations between the havon VEVs involving different 
A couplings which leads to an additional constraint on the phases of the respective VEVs. 
Exact values of these A are not specihed; it sufhces that they are not equal. 


Furthermore, the primary role of the SU(5) adjoint helds Z 2 and which couple to 
is in the GUT breaking mechanism (as discussed in Section 4.1). Their phases are hxed 
separately by other superpotential terms. 


We begin the analysis of these terms by noting they do not ahect the alignments of 
the havons 0. The corresponding P-terms for each held Pij produces a set of coupled 
equations that admit a solution where none of the A, O, and P helds but all the havons 
obtain a VEV. Omitting the (real, positive, 0{1)) A coefficients, these VEVs have the 
structure: 


Ve 

~ iff 

Vf, 



Vr 

~ {VZ2 

Vl 

~ iyz^vzsY 

2 vl 

V2 

~ {vz 2 Vz 3 )~^M~i vl 

V 4 

~ {VZ2VZ3)~ 

6 Ms U3 

V 5 

~ M~^ U3 

Ve 

~ {VZ2VZ3)~ 

1 , .r 2 9 

eM 3 vl 

Vdi 

~ (uZjUZgjsiffS Ug 

ve2 

~ (VZ2UZ3)" 

1,^11 

6 M 3 Vg 

Vo 

= vp = va = 0. 
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Regarding the magnitudes of the VEVs, two comments are in order. We assumed above 
that M sets the scale of the VEV of (/)e, which is in contradiction with our previous 
assumption that it be (9(10“^)M. This violates our simplifying assumption that all mass 
scales are equal, and demonstrates that some spectrum of mass scales is in fact required 
in this model. As for the VEV U 3 , as discussed earlier, it is driven to a specihc scale A 3 
radiatively 25 . Writing pi = arg(uj), this VEV structure gives (up to multiples of vr) the 


phase relation 


27rn 1 
Pi = - 6 


(3.5) 


where n is an integer, and similar relations for the other flavons as linear combinations 
of p 3 , {pz 2 + PZi) and multiples of 27r/3. This is an important equation since it Exes the 
relative phase ps — p^ = Patm — Psoi in terms of |(pz 2 + Pzs)- As discussed in Section 
PZ 2 + PZ 3 = where we also establish that p^ = for integers k, k'. From Eq. 


4.1 


2.16 


p = — pg + 2 (patm— Psoi), SO we couclude that p is one of the nine complex roots of unity. 


GUT scale symmetry breaking, proton decay and 
the strong CP problem 


In this section we discuss the aspects of the model related to grand unihcation, starting 
with how the R-symmetry and the GUT gauge group are spontaneously broken. We refer 
to the superfields involved as the scalar sector; they are listed in Tables and We then 
describe the details of the MP mechanism, and hnish this section with an analysis that 
justihes the absence of dangerous proton decay operators in the model. 


4.1 SU{b) and Zf breaking 

As previously discussed, the T messengers form pairs; their mass scale, unprotected by 
any symmetry, is near the highest scale of the theory, which we represent generically as 
M. The GUT breaking superpotential with non-renormalisable terms is theiUI 

IUguT = Zi ^MA24 + + Z 2 ^ 

+ ^3 (^5-^24 + ^6^3) ■ 

We have hve GUT adjoint superhelds, three of which (the Z helds) are charged by 2 and 
two (A 24 and II 24 ) by 0 under the R-symmetry. Also appearing in Wgut is the Majoron 

renormalisable term of the form Z 2 i? 24 n 2 , allowed by the symmetries, has been dropped to make 
the discussion more transparent. This term mixes the VEVs of the GUT breaking scalars with the ones 
in the MP mechanism so they should be naturally around the same scale (M ^Mqut)- Beyond this, its 
practical effect is minimal: the fields obtain VEVs with or without this term. Since the VEVs get very 
complicated when this “mixing” term is included, we ignore it for simplicity, simply bearing in mind 
that VEVs from both sets of fields are related. 
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the GUT singlet field which we have seen is involved in giving mass to several SM 
fermions and whose VEV breaks lepton number by giving the right-handed neutrinos 
their Majorana mass. The supressions of the non-renormalisable terms in Eq. 4T come 
precisely from the mass of the T messengers, as displayed in Fig. 


I 


I 


(a) 


H-21 


T« T? 

(b) 


Figure 4: Diagrams for non-renormalisable potential terms. 


VEgut has a non-trivial minimum: 

3V2A2A3Af 
t2V3AiA^/^ 
SUSA^/^AsAg/" 
tAiA^" 

SVaA^/^^A.Aj/^Ay® 


VZr = 


vz. = 


vz. = 


M 


M 


M 


VH24 = 


^A24 


( 0 ^ = 


AiA 4 ^"A 


1/3X1/3 


6 


-M 


A^^AsAs 


2/3 

A 2 AIA 5 

2^,,3 
Ay^A, ’ 


M 


(4.2) 


where all the adjoint scalars get a VEV of the form ($ 24 ) = 'f 4>24 diag(2, 2, 2, —3, —3). By 
themselves, the F-terms associated with PUgut also allow a trivial minimum where the 
magnitude of each VEV vanishes. But after SUSY is broken and we consider the effects of 
the small contribution from radiative breaking 23 to the scalar components of the GUT 


breaking superfields (as we did in Section]^ for the ^44 breaking fiavons), the stationary 
point with vanishing magnitudes is no longer a minimum due to the radiatively induced 
negative squared mass term. To a very good approximation the true minima are given 


by the magnitudes in Eq. |4.2[ which are now a lower energy state than the trivial F-term 
solution. 

We conclude that Eq. can generate GUT and R-symmetry breaking at high scale, 
with Zf broken to (standard R-parity preserved) by the Zi VEVs. Because Zf is 


broken at a high scale, the no-go theorem from [ 2 ^ does not apply to our model and we 
verified that all the components of the SU{5) adjoints acquire GUT scale masses. 

A slightly unappealing issue with Wgut is that the minimum requires some non-(!2(l) 
choice of A parameters if we are to obtain a hierarchy between the VEVs of F 24 and 


A 24 , and an appropriate value for {^/M as shown in Eq. 2.22 These requirements come 
from the successful fit to up and down quark and charged lepton masses, as discussed 


in Sections 24^, 2 ^ and 2 G, and partly also for the /i term, as will be discussed shortly. 
However, since the messengers will in general have different masses (recall we set them all 


equal to M only for simplicity), the A parameters need not be as hierarchical as Eq. 4.2 
appears to indicate. For example, if the masses of messengers E are slighly larger than 
the GUT scale masses of messengers T, this would allow all A to be 0{1). 
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We note also that, although we are considering a situation where the superpotential 
parameters (the M and A couplings) are real due to CP conservation, the VEVs of the 
GUT breaking scalars may be complex, since they depend on order roots of real 
numbers. As has been shown in Sections 2 ^ and 3^, the phases of the helds Z 2 and 
Z 3 are relevant for establishing the physical phase rj in the neutrino mass matrix, which 
controls neutrino masses and mixing. We see immediately that for integer k, 

i.e. one of nine roots of unity. While pz 2 and pz^ individually can be any of six roots, 
originating in the factor their product ^’ 2^3 cancels this factor such that the largest 
root is a third, giving pz 2 + Pza = for integer k'. 


4.2 Doublet-triplet splitting, Higgs mixing and the /i term 


Given that we have a number of GUT representations containing (SM gauge group) SU{2) 
doublets and triplets we turn now to a brief discussion of how doublet-triplet splitting 
is achieved in this model. Although one could alternatively introduce the double MP 
mechanism 27 as demonstrated in (T^, here we limit ourselves to describing how the 


MP mechanism is implemented with the helds listed in Table 
We have a superpotential 


Wn = ni 


^A7n^ -1- mt\2 t 



(4.3) 


which gives If, the 75s of 57/(5), their VEVs 






vn2 



(4.4) 


which are aligned with the SM singlet inside the SU{5) 75. The non-renormalisable term 
in Wn comes from the diagram in Fig. 


I 

I 

I 

I 




I 

I 

I 


' nj 
I 
I 
I 


n. 





Ho 


Figure 5: Diagram for non-renormalisable term in the 75 potential. 


With Eq. |4.4| , the MP mechanism proceeds through the superpotential: 

IUmp ~ H-^n,U2 + H5n2lli + 

-|- 7745723112 -|- MH 


(4.5) 
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The very high order non-renormalisable term at the end arises through the S messengers 
already listed in Ta ble (half of which participate in the up quark Yukawa terms, as 
discussed in Section 2.1 and illustrated in Fig. [^. Strictly speaking this term does not 
participate in splitting the masses of SU{2) doublets and triplets, rather it is the source 
of the /i term in our model, as shown below. 


The terms in ITmp generate mixing between the 45s and 5s of SU{5). The mass matrix 
for the SU{2) triplets contained in the 45s and 5s is: 



( 0 


vn2 

0 \ 

—T 

vu,e 

M 

0 

vn2 


vni 

0 

(0 

0 


\ 0 

0 

0 

M / 


(4.6) 


3^={3{H-,) 3{H,-,) 3 (^ 2 ) 3 ( 124 )) 

3^ = {3iH,) 3iH,,) 3(12i) 3(123)), 

where once again ^ = v^/M. Taking (ni^ 2 ) ~ M, the eigenvalues of this mass matrix are 
all of order M (i.e. at the GUT scale), leading us to conclude there are no light SU{2) 
triplets. Conversely, for the doublets we have the matrix: 


(mi) 



t-nA (2{H,) \ 
M ) ' 


(4.7) 


It is clear that were it not for the determinant of this mass matrix wonld vanish. 
We may rotate to the basis of the MSSM Higgs doublets s and a pair of very heavy 
doublets 


(2{Hl)\ 

\ms)) 

{mi)\ 


1 / 2(1 1 ) 





H.. 


(4.8) 


The usual MSSM term comes from this mechanism with: 


M 


2.22 


(4.9) 


we see that 


where vpi provides the necessary breaking. Using the £t from Eq. 

~ 1.6 X 10“^°Mgut- If we choose the couplings at the vertices of the tower that 
generates the term to be ~ 0.5 we may get a term /i ~ 0(10^ — 10^) GeV without any 
hne-tuning. 


4.3 Proton decay 

A classic problem in GUT theories, and in particular those based on SU{5), is the pre¬ 
diction of excessively fast proton decay. The most dangerous processes come from the 


21 







“dimension 5” operators TTTF (for a discussion of dimension 6 operators we refer the 
reader to [^). The “dimension 5” operators are forbidden by the symmetries of the 
model, but related higher order operators of the following form are allowed: 


T^T^nF 


Zcj) 

W 



(4.10) 


where the extra superhelds shown are needed for such terms to be invariant under the 
symmetries. Since we are working with the renormalisable theory, in order for this type 
of effective term to be present at the GUT scale, there must be messengers allowing 
them. In this case, an analysis of the SU{5) index structure revels there should either be 
messengers that are SU{5) 10, or SU{5) 5 that are also charged under Zf. As one can 
conhrm from Table our model has neither: 10 messengers were not used, and the 5 
messengers are all neutral under Zf. We conclude therefore that our symmetry content, 
together with the existing set of messengers, do not allow any such GUT scale suppressed 
operators that would lead to excessively fast proton decay to be generated. The operators 
in Eq. 4.10 may in principle be generated by physics at the Planck scale, with the scale 
M replaced by the Planck mass, leading to highly suppressed proton decay. 


4.4 Strong CP problem and the Nelson-Barr resolution 


We hrst recall the strong GP problem, namely that the physical angle 6 = 6 *qcd — 
where 6 'qcd multiplies the topological gluon term {g^/327i‘^)GG and O g = arg detfU^U'^), 

The origin 
It is inter- 


20 


,|28 


is limited to be 6* < 10“^*^ by the non-observation of the neutron EDM 
of such a small number, 6 < 10“^°, is commonly called the strong GP problem 
esting to compare this with the GP violation related to the weak interaction in the quark 
sector; the relevant quantity is the Jarlskog invariant det[y“yU,yrfydt], ^pich, 

when compared to data, is required to be non-vanishing, and indeed in the standard 
parameterisation, requires a large phase angle 5'^ ~ 1. 

It turns out that our model resolves the strong GP problem without relying on the 
introduction of axions. Unlike the axion solution, which requires a global U{1) symmetry 
with a colour anomaly, we shall rely on the fact that the high energy theory conserves 
GP, ensuring that ^qcd = 0. GP is then spontaneously broken in such a way as to yield 
(59 ~ 1 while maintaining 6 < 10“^° and in particular 6q < How it achieves this 

feat can be seen from Eq. 2.18 where is real, while the structure of in Eq. 2.19 gives 
it a real determinant. This is due to the lack of a Yukawa term Y 21 , meaning the coupling 
¥^2 (which is the only complex Yukawa coupling) does not appear in the determinant of 
Y'^. Therefore there are no contributions to 6q even after spontaneous GP breaking. This 
is similar to the Nelson-Barr mechanism 29,30 , where the triangular form of Yukawa 


matrices was proposed, although in our model 9q vanishes due to the triangular form of 
only, with Y“ being non-triangular and real. 


For a successful resolution of the strong GP problem, we must ensure that no higher order 
corrections to the Yukawa matrices arise which would violate the bound 6q < 10“^°. The 
main focus is on the Yukawa coupling Y 21 which is zero at leading order but which may 
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in principle receive higher order corrections, violating the triangular structure. However 
in our model such higher order corrections are absent at the held theory level with the 
specihed messenger sector. This entry in the Yukawa matrix would arise from the coupling 
of the bilinears T 2 H^ 45 to the bilinear (l)eF. Since these terms are non-renormalisable, we 
require messengers to form them. The messengers that could produce such terms are the 
Xi helds in Table With these messengers, the only allowed connection to is 
(contributing to Yf^), thus forbidding the Y 2 I term, even when allowing for all higher- 
order corrections. Therefore the specihed model has no strong CP violation arising from 
Y 21 since the required operators are not generated at the held theory level. 


It is also important to consider the ehect of higher order corrections arising from the 
Planck scale, since such operators only have to respect the symmetries of the model, and 
do not require the specihed messenger sector to generate them. The biggest contribution 
would come from the termEE 

(4.11) 


‘Mp 


With a general choice of phase, such a term would lead to 9g ~ 10“"^ which is far too big. 
However the contribution to 6 q may be avoided by a judicious choice of GUT breaking 
phases. As stated in section 2.4, the physical phase in the down quark Yukawa matrix 
is C = P? — ‘^Ph 2 a ~ PA 24 - The new Planck suppressed term has a phase C,' = —p^ + 
2pA24- Choosing a relation between phases 2 ph2a = PA 241 then C, = —C' and this way 
the contribution to 9q vanishes. This happens for one in three cases. The next biggest 
contribution comes from a term 

^ ^ (4.12) 




45 


Ml 


giving 9q ~ which is several orders of magnitude below the current experimen¬ 

tal bound. Any other Planck suppressed terms allowed by the symmetries are further 
suppressed so we need not consider them. Therefore the model may resolve the strong 
CP problem even in the presence of Planck scale operators controlled only by symmetry. 
Finally, extra contributions may come from SUSY breaking terms. If we assume that 
there is no extra CP violation in this sector, which is controlled by the spontaneously CP 
violating flavons, the SUSY flavour problem is under control and such contributions to 9 


are also expected to be negligible 31 


5 The leptogenesis link 


The link between leptogenesis and the PMNS matrix was hrst studied for sequential 
dominance in 
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In the original form of CSD, the columns of the Dirac mass matrix 
in the flavour basis were orthogonal to each other and consequently the CP asymmetries 
for cosmological leptogenesis f 


vanished 34-37 . In this model, leptogenesis does not 


vanish since the columns of the Dirac mass matrix in the flavour basis are not orthogonal. 
^^This term would also give a contribution to lepton angles of 0(10“^) which is negligible. 
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Interestingly, since the seesaw mechanism in this model with two right-handed neutrinos 
only involves a single phase t] = both the leptogenesis asymmetries and the neu¬ 

trino oscillation phase must necessarily originate from this phase, providing a direct link 
between the two CP violating phenomena in this model. 


Following the arguments in (34], the produced baryon asymmetry Yb from leptogenesis 


in the seesaw model in Eq. 2.11 satishes 


Yb oc ± sin 77, 


(5.1) 


where the “-I-” sign applies to the case Matm ^ Mioi and the ” sign holds for the case 
Msoi dPatm- Since the observed baryon asymmetry Yb is positive, it follows that, for 
d^atm ^ dPsoh we must have sin 77 to be positive, while for Mgoi Matm we must have 
sin ?7 to be negative. We have seen that positive rj is associated with negative 6 ^ and 
vice versa. Although the global hts do not distinguish the sign of t], the present hint 
that (5^ ~ — 7 r /2 would require positive t], then in order to achieve positive Yb we require 
Matm "C Msoi, which is natural in our model, corresponding to “light sequential domi¬ 
nance” , where successful leptogenesis may be achieved in the two right-handed neutrino 
model as discussed in 
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6 Conclusion 


We have presented here a fairly complete realisation of an SU{5) GUT flavoured with 
A 4 , which leads to the MSSM plus two right-handed neutrinos below the GUT scale. 
The A 4 family symmetry unihes the three families of 5-plets F and its vacuum alignment 
determines the Yukawa matrices. In addition a Zg x Zg symmetry provides the mass 
hierarchy and controls spontaneous CP violation in both the quark and lepton sectors 
while a Z^ symmetry is broken to Z^, identihed as the usual R-parity. Proton decay 
is under control in this model, with the symmetries forbidding dangerous dimension-5 
operators, and similar (but higher order) operators being very suppressed. The strong CP 
problem is resolved in a similar way to the Nelson-Barr mechanism. The model is highly 
predictive and satishes many distinct (and non-trivial) phenomenological requirements. 

Imposing CP at the high scale is an important feature of the model. If we do not impose 
CP then all couplings become complex, leading to all VEVs having undetermined phases. 
In particular the phase rj, present in both neutrino mixing and leptogenesis, is no longer 
restricted to a discrete choice. However the link between leptogenesis and low energy 
phenomenology remains. On the other hand we would no longer solve the strong CP 
problem. 

We highlight the ubiquitous nature of the havon held ^ across all the sectors of the model: 
it triggers spontaneous CP violation in both the quark and lepton sectors, generates up- 
type quark mass hierarchies and CKM mixing, explains the smallness of down quark 
and electron masses and breaks lepton number, providing the hierarchy between solar 
and atmospheric right-handed neutrino masses. In addition, the ^ held is responsible for 
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generating the small (complex) fi term. The phase of the ^ field VEV also contributes 
to the relative phase 77 appearing in the neutrino mass matrix. This phase, the only one 
appearing in the neutrino mass matrix and in the formula for the baryon asymmetry of 
the universe, provides a direct link between the PMNS phase 5^ and leptogenesis. 

We emphasise that the entire PMNS matrix is predicted without any free parameters, 
up to nine choices for 77 , where we select 77 = +27r/3, since it is preferred by comparing 
CSD3 to data. The required vacuum alignments are provided from ^ 4 . The model 
effectively serves to yield the CSD3 scheme (with two right-handed neutrinos) within a 
fully working and viable SUSY GUT of flavour in which all quark and lepton masses and 
mixings are successfully described. Indeed the model provides an excellent £t (better 
than one sigma) to all quark and lepton (including neutrino) masses, mixing and CP 
violation. All fermion mass hierarchies are understood in the sense that purportedly 0{1) 
couplings indeed contain no strong hierarchies. However the most immediate predictions 
of the model are those of CSD3 with two right-handed neutrinos and 77 = -|-27r/3, namely 
a normal neutrino mass hierarchy with mi = 0, a reactor angle of 6[^ ^ 8.7°, a solar 
angle 0(2 ~ 34°, close to maximal atmospheric mixing ^23 ~ 46° and almost maximal 
leptonic CP violation, with an oscillation phase ~ —87° consistent with the current 
hint S’- ~ —7r/2. 

The reason why the held content is so large is that the model is fairly complete. In 
particular it is renormalisable at the GUT scale, which requires a large explicit held 
content including many heavy messenger superfields. It also addresses many aspects 
relevant both to a GUT and to family symmetry models (stopping short of discussing the 
details of SUSY breaking and its string theory completion). In particular, the A 4 , SU{5) 
and i?-symmetry symmetry breaking sectors all require large held content. For example, 
the GUT symmetry is broken by an explicit superpotential at the GUT scale, including 
doublet-triplet splitting via a missing partner mechanism (leaving no light exotic degrees 
of freedom at the low scale), Higgs mixing and the origin of the MSSM /i term of the 
right order of magnitude. 

Despite the many successes of the model, there are inevitably several important issues 
that he beyond the scope of this paper. To take one example, we do not discuss GUT scale 
threshold corrections, which will be important in maintaining successful gauge coupling 
unihcation in the presence of many fields, including colour triplets, at the GUT scale. 
In fact all the additional superfields in non-trivial representations of the gauge group 
may have masses at or above the GUT scale. Another important issue is that of the 
low energy superpartner spectrum in this model. While we expect SUSY induced flavour 
changing to be under control for the 5-plets, which are unified into an A 4 triplet, this is 
not the case for the 10 -plets Tj which are singlets of A4, leading to flavour violation in 
the super-GKM basis. It would be interesting to study this in the future. 
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A A4 basis convention 


In the basis we are using (see 17 for more details), one has the following Clebsch-Gordan 
rules for the multiplication of two triplets, 3 x 3 = 1 -|- 1' -|- 1" -|- 3i -|- 82 , 


(ab)i = aibi 02^2 + 03^3 
(ab)p = aibi + coa2b2 + co^asbs 

(ci6)i" = (iibi -|- ca^ci2^2 T ^®3^3 (-^'1) 

(06)31 = (0263,0361,0162) 

(06)32 = (0362,0163,0261) 


where the components of the two triplets are given by o = (01,02, 03) and 6 = (61, 62, 63), 
and = 1. 
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